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TECHNICAL  RELEASE 
NUMBER  15 

COMPUTATION  OF  WATER  SURFACE  PROFILE 
AND  RELATED  PARAMETERS  BY  ES-CHARTS 


CHAPTER  1 - DESCRIPTION  OF  ES-CHARTS 

The  theory,  definitions,  and  nomenclature  used  in  this  Technical  Release 
are  presented  in  the  National  Engineering  Handbook,  Section  5,  Hydraulics. 

The  procedure  for  water  surface  profile  determination  is  a modification 
of  Escoffier’s1  method  for  natural  channels.  The  modification  includes 
the  effect  of  changes  in  velocity  head  on  the  profile,  and  the  charts  are 
a slightly  different  plot  than  used  in  Escoffier’s  original  presentation. 

The  procedure  presented  is  useful  in  the  design  of  prismatic  channels. 

It  readily  shows  the  effect  on  the  profile  of  changes  in  the  tailwater 
condition,  bottom  slope,  Manning’s  roughness  coefficient,  and  bottom  width 
of  the  channel.  These  charts  are  useful  in  determining  the  profiles  asso- 
ciated with  various  steady  discharges.  They  are  not  applicable  for  un- 
steady flow  conditions. 

The  charts  are  prepared  for  rectangular,  trapezoidal,  and  circular  pris- 
matic channels.  A catalog  of  the  various  charts  available  is  given  on 
page  1-2.  Generally,  six  sheets  are  used  to  present  the  parameters  for 
each  bottom  width  of  a channel  with  a particular  side  slope  z:l.  A lap 
of  about  one -third  of  the  depth  increment  is  shown  from  one  sheet  to  the 
next.  The  sheets  are  so  arranged  that  the  curves  on  one  side  of  any 
sheet  do  not  lap  the  curves  on  its  reverse  side. 

Numerous  combinations  of  grid  sizes  were  used  to  permit  a reasonable  de- 
gree of  accuracy  of  the  graphical  solution.  The  values  given  by  the 
charts  were  obtained  entirely  by  the  use  of  an  electronic  digital  com- 
puter. Data  were  computed  at  depth  intervals  equal  to  or  smaller  than 
the  scale  graduations  of  the  chart.  The  original  drawings  of  charts 
were  twice  the  linear  size  of  the  multilithed  prints. 

Each  sheet  has  two  independent  line  charts  (see  page  1-3)-  The  vertical 
line  chart  on  the  left  side  of  the  sheet  is  composed  of  the  two  scales 
Qc,d  and  d*  The  d-scale  is  the  same  as  the  scale  on  the  axis  of  ordinates 
of  the  central  chart . 


1 


Escoffier,  Francis  F.,  Graphical  Calculation  of  Backwater  Eliminates 
Solution  by  Trial,  Engineering  News -Record,  June  27,  19^-6. 
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CATALOG  OF  ES -CHARTS 

FOR  COMPUTATION  OF  WATER  SURFACE  PROFILES 


z 

* 0:1 

* 1:1 
* 1^:1 
* 2:1 

* 2j:l 

* 3:i 
***  4:1 
***  6:1 

**  circular 


b 

Sheet 

Numbers 

Depth  Range 

b 

Sheet 

Numbers 

Depth  Range 

1 

1 

- 

6 

0.2 

- 10.0 

70 

163 

- 

168 

0.2 

- 30.0 

2 

7 

— 

12 

0.2 

- 10.0 

75 

I69 

— 

174 

0.2 

- 30.0 

3 

13 

— 

18 

0.2 

- 10.0 

80 

175 

— 

180 

0.2 

- 30.0 

4 

19 

— 

24 

0.2 

- 10.0 

85 

181 

- 

186 

0.2 

- 30.0 

5 

25 

— 

30 

0.2 

- 10.0 

90 

187 

— 

192 

0.2 

- 30.0 

6 

31 

— 

36 

0.2 

- 15.0 

95 

193 

— 

198 

0.2 

- 30.0 

7 

37 

— 

42 

0.2 

- 15.0 

100 

199 

— 

204 

0.2 

- 30.0 

8 

43 

— 

48 

0.2 

- 15.0 

110 

205 

- 

210 

0.2 

- 30.0 

9 

49 

- 

54 

0.2 

- 15.0 

120 

211 

- 

216 

0.2 

- 30.0 

10 

55 

- 

60 

0.2 

- 15.0 

130 

217 

— 

222 

0.2 

- 30.0 

12 

61 

- 

66 

0.2 

- 20.0 

140 

223 

- 

228 

0.2 

- 30.0 

14 

67 

- 

72 

0.2 

- 20.0 

150 

229 

- 

234 

0.2 

- 30.0 

16 

73 

- 

78 

0.2 

- 20.0 

160 

235 

— 

240 

0.2 

- 30.0 

18 

79 

- 

84 

0.2 

- 20.0 

170 

241 

- 

246 

0.2 

- 30.0 

20 

85 

- 

90 

0.2 

- 20.0 

180 

247 

- 

252 

0.2 

- 30.0 

22 

91 

— 

96 

0.2 

- 25.0 

190 

253 

— 

258 

0.2 

- 30.0 

24 

97 

— 

102 

0.2 

-25.0 

200 

259 

— 

264 

0.2 

- 30.0 

26 

103 

- 

108 

0.2 

- 25.O 

210 

265 

— 

270 

0.2 

- 30.0 

28 

109 

— 

114 

0.2 

- 25.O 

220 

271 

— 

276 

0.2 

- 30.0 

30 

115 

- 

120 

0.2 

-25.O 

230 

277 

— 

282 

0.2 

- 30.0 

35 

121 

- 

126 

0.2 

- 30.0 

240 

283 

— 

288 

0.2 

- 30.0 

40 

127 

- 

132 

0.2 

- 30.0 

250 

289 

— 

294 

0.2 

- 30.0 

45 

133 

— 

138 

0.2 

- 30.0 

260 

295 

— 

300 

0.2 

- 30.0 

50 

139 

— 

144 

0.2 

- 30.0 

270 

301 

— 

306 

0.2 

- 30.0 

55 

145 

— 

150 

0.2 

- 30.0 

280 

307 

— 

312 

0.2 

- 30.0 

60 

151 

- 

136 

0.2 

- 30.0 

290 

313 

— 

318 

0.2 

- 30.0 

65 

157 

— 

162 

0.2 

- 30.0 

300 

319 

— 

324 

0.2 

- 30.0 

* Available  for  b ^ 75* ** ***  Available  for  80  ^ b ^ 300  by  special  request 
Delivery  of  special  request  orders  is  usually  slower  than  available 
stock  material. 

**  Available  after  January  1962 

***  Not  to  be  drawn  by  Design  Section.  A limited  supply  of  computations 
(1  ^ b £ 300)  is  available  for  plotting  charts. 
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The  vertical  line  chart  located  on  the  right-hand  side  of  the  sheet  is 

”nQn,  d" 

composed  of  the  two  scales  


lSo1'! 


and  d.  This  d-scale  is  also  the  same 


as  the  scale  given  on  the  axis  of  ordinates  of  the  central  chart. 


These  line  charts  may  be  used  to  determine  various  hydraulic  parameters. 
The  determination  of  some  of  these  parameters  will  require  a mental  re- 
designation  of  the  scales  of  the  line  charts.  This  will  not  be  diffi- 
cult if  the  relations  of  the  various  parameters  axe  well  understood. 

Such  redesignations  are  given  under  Chapter  2 "Determination  of  Hydraulic 
Parameters" . 

The  use  of  the  central  chart  and  the  two  orthogonal  lines  whose  scales 

0 

have  the  labeled  values  m and  — will  be  discussed  later. 

m 


CHAPTER  2 - DETERMINATION  OF  HYDRAULIC  PARAMETERS  FROM  LINE  CHARTS 


The  procedures  for  determining  various  hydraulic  parameters  are  given 
in  the  following  numbered  paragraphs.  The  examples  1-10  on  pages  ^-1 
through  A -5  are  numbered  to  correspond  with  the  numbered  procedures 
given  here. 

Procedure  1. 


The  critical  discharge  Qc?£  of  a given  channel  corresponding 
to  a depth  d is  read  directly  from  the  left-hand  line  chart. 

Procedure  2. 


The  critical  depth  dc>Q  of  a given  channel  corresponding  to  a 
discharge  Q is  read  directly  from  the  left-hand  line  chart. 
When  the  discharge  Q is  given,  the  parameters  of  this  line 
chart  are  mentally  redesignated  as  Q and  dCjQ. 

Procedure  3- 


The  normal  discharge  Qnt(j  of  a given  channel  corresponding  to 
a depth  d is  determined  by  the  use  of  the  right-hand  line  chart. 

a.  For  the  depth  d read  the  corresponding  value  of 

~nSn,al  . 


s 1/2 
lPo  j 


b . Compute  the  value  of 


n 


s x/2 

Lro  J 


for  the  given  channel,  or 


find  the  value  of 


n 


lPo17' j 


from  ES-138. 


Divide  the  value  of 


~nQn,  d 


s V2 

LPO  -J 


read  from  the  line  chart 


by  the  value  obtained  for 


n 


So1/2 


Procedure  4. 


The  normal  depth  dn^  q of  a given  channel  corresponding  to  a 

discharge  Q is  determined  from  the  right-hand  line  chart.  When 
the  value  of  Q is  given  and  the  value  of  dn  q is  to  be  deter- 
mined, the  parameters  of  this  chart  are  mentally  redesignated 

" nQ  1 n , 
as  and  dn^  q. 


1/2 


LrO  J 


a.  Compute  the  value  of 


n 


Lso1^ 


for  the  given  channel 


2-2 


t 


or  find  the  value  of 


n 


s x/2 
bo 


from  ES-1^8. 


nQ 

n 

■cn 

O 

H 

lJL 

_So1/2_ 

h . Compute 

c . For  the  computed  value  of 
ponding  normal  depth  dn^  q . 


Q. 


nQ 


q x/2 
LrO  J 


read  the  corres- 


Procedure 

The  normal  slope  sDi of  a given  channel  corresponding  to  the 
discharge  Q and  depth  d is  also  determinable  from  the  right- 
hand  line  chart.  When  the  values  of  Q and  d are  given  and  the 
corresponding  normal  slope  sn  is  to  be  determined,  the  par- 
ameters of  the  line  chart  are  mentally  redesignated  as 

nQ 


and  d. 

_(sn,Qd)1/;: 

a.  For  the  depth  d read  the  corresponding  value  of 
nQ  ~j 

lsn,Qd)1/2 

b.  Method  1.  Compute  nQ,  and  divide  by  the  value  of 


nQ 


isn,Qd)1/! 
That  is. 


Square  the  result  to  obtain  sn^Qd< 


3n,  Qd 


nQ 

nQ 

>n,  Qd ) 1 /2 

Method  2.  Divide  the  value  of 


nQ 


by  Q to 


determine 


n 


L(sn,Qd)17  J 
parameters  [ sD,  n,  and 


n 


L<sn,Qd)1/2_ 
Mentally  redesignate  the 


LSo1/2j 


jof  ES-138  as 


sn,Qd>  n>  and 


n 


iQn.QdfQj 


For  the  given  n and  determined 
sn, Qd . 


n 


isn,Qd)1;l 


, read 


Procedure  6. 


The  critical  slope  sCi£  of  a given  channel  corresponding 
to  the  depth  d is  determined  by  the  use  of  both  the  right- 
and  left-hand  line  charts.  The  parameters  of  the  right - 

nQr 

hand  line  chart  are  mentally  redesignated  as 


*c,  d 


_(sc,d)1/2_ 


and  d when  the  critical  slope  corresponding  to  a depth  is 
to  be  determined. 


a.  Read  the  values  of  Qc^  and 


nQ 


ClA. 


correspond- 


.(Bc,d)1/S 

ing  to  the  depth  d from  the  two  line  charts. 
Method  1.  Compute  nQc^d,  and  divide  by  the  value 


of 


nQ, 


SjA. 


_(sC,d)1/2_ 

That  is 


Square  this  result  to  obtain  sc^. 


;c,  d 


nQc,a 

r n^c,d  1 

1 

>c,d)1/2_ 

2 


Method  2.  Divide  the  value  of 


n 


nQ 


c,  d 


Qc.d  to 

L(sc,d)1/2. 

. Mentally  redesig- 


n 


3o 


/2 


of  ES-1J8  as 


obtain  the  value  of 

.(sC,d)1/2 

nate  the  parameters  ^s0j  n,  and 
(sC,a>  n,  and  

V L(sc,d)x/2Jy 

sc,d  from  ES-I38.  For  the  given  n and  determined 


Read  the  value  of 


n 


(Bed)1/2. 


y read  sc^. 


Procedure  7« 


The  critical  slope  sCjq  of  a given  channel  corresponding  to 

the  discharge  Q is  determined  by  the  use  of  both  the  right - 
and  left-hand  line  charts.  The  parameters  of  the  left-hand 
line  chart  are  mentally  redesignated  as  Q and  d-c,Q>  an(^  the 


2-4 


parameters  of  the  right-hand  line  chart  are  mentally  redesig- 

and  dn^Q. 


a.  For  the  discharge  Q read  the  corresponding  value  of 
dCjQ.  By  the  definition  of  critical  slope  observe 
that  dc^Q  = dn>Q. 


nated  as 


nQ 


(sc,q) 


-/2 


b. 


c . 


For  the  value  of  dc^Q 
nQ 


value  of 


_( sc,  q)  1/2_ 


= dn^Q,  read  the  corresponding 
from  the  right-hand  line  chart. 


Method  1. 
nQ 


Compute  nQ,  and  divide  by  the  value  of 
, and  square  the  result  to  obtain  the 


_(sc,  q)  1/2_ 
critical  slope. 

That  is 


;c,Q  ~ 


nQ 

nQ 

_(sC,q)1/2_ 

Method  2.  Divide  the  values  of 


nQ 


_( sc,  q)  1/2_ 


by  Q 


to  determine 


n 


_(sc,q)1/2 

Mentally  redesignate  the  parameters  (s0,  n,  and 


n 

s 1/2 
LrO 


of  ES-I38  as-  n,  and 

given  n and  determined 


n 


lKq)17! 


For  the 


n 


_( sc,  q) 


/2 


, read  sc^q. 


Procedure  8. 

dhf 

The  rate  of  friction  loss  s = -j-j—  corresponding  to  the  dis- 

charge  Q and  depth  d is  determined  by  the  use  of  the  right- 
hand  line  chart.  Theory  shows  this  rate  is 


s = 


dhj 

dT 


= Sr 


Q 


-1 2 


_9n,  d 


where  s is  dimensionless  or  has  the  units  foot  pound  per  foot 


2-5 


length  of  channel  per  pound  of  water. 

a.  For  the  depth  d read  the  corresponding  value  of 
~nQn,d~ 


LSo1/2_ 
t> . Divide 


<d 


~nQn,d" 

;o 


s,1/2 


by  n,  and  square  the  result  to  obtain 


so 

Divide  Q2  by 


Q§,d" 


_ so 


to  obtain  s = sf 


P-’ 

P,  d 


Procedure  9» 


The  direction  of  the  solution  line  is  required  before  water 
surface  profiles  can  be  determined.  It  is  also  used  to  deter- 
mine velocity  head,  specific  energy  head,  and  the  total  friction 
loss  between  the  two  sections  used  to  determine  a depth  of  flow. 
The  direction  of  the  solution  line,  for  a given  sheet,  depends 
only  on  the  discharge  Q. 

The  direction  of  the  solution  line  is  determined  from  the  two 
orthogonal  lines  found  along  the  bottom  and  extreme  left-hand 
side  of  the  sheet.  Each  of  these  orthogonal  lines  have  a 
single  scale.  (See  page  1-5*)  The  value  of  m selected  must  be 

Q 

used  in  determining  the  value  of  — . On  any  one  sheet  and  for 

a given  value  of  Q,  the  direction  of  the  solution  line  is  the 
same  for  any  selected  value  of  m. 


The  following  procedure  is  used  to  determine  the  direction 
of  the  solution  line. 

a.  Select  the  smallest  value  of  m for  which  the  corres- 

Q Q 

ponding  value  of  — can  be  read  on  the  — scale. 

m m 

b.  The  direction  of  the  solution  line  corresponding  to 
the  discharge  Q is  the  same  as  that  vector  drawn  from 

Q 

the  point  m to  the  point 

The  tangent  of  the  angle  which  the  solution  line  makes  with 


the  axis  of  abscissas  is 


_2g 


The  direction  of  the  solution 


line  corresponding  to  the  discharge  Qc^  is  parallel  to  the 

tangent  of  the  curve  at  the  depth  d. 
a 


') 
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Procedure  10. 

The  converse  situation  of  determining  the  discharge  Q,  corres- 
ponding to  a given  direction  of  the  solution  line  frequently 
arises . 


a.  Select  the  smallest  value  of  m on  the  m-scale  such 
that  a line  having  the  given  direction  will  intersect 

Q 

the  — scale  when  drawn  from  the  point  m. 
m 

Q 

b . Read  the  value  of  — at  the  intersection,  and  multiply 

by  the  selected  value  of  m to  obtain  the  correspond- 
ing value  of  discharge  Q. 


3-1 


# 


CHAPTER  3 - DETERMINATION  OF  SPECIFIC  ENERGY  HEAD, 

VELOCITY  HEAD,  AND  WATER  SURFACE  PROFILES. 

The  method  of  determining  various  parameters  and-  the  water  surface  pro- 
files by  use  of  the  central  chart  are  given  by  the  following  numbered 
procedures.  The  examples  on  pages  4-5  through  4-17  illustrating  the 
use  of  the  procedures  are  numbered  to  correspond  with  the  numbered 
procedures  given  here. 

Procedure  11. 


4 


The  specific  energy  head  is  the  vertical  distance  between 
the  energy  grade  line  and  the  bottom  of  the  channel.  Its 
value  is  often  used  in  hydraulic  calculations.  The  graph- 
ical solution  for  the  specific  energy  head  H^q^  of  a 

given  channel  corresponding  to  the  depth  d and  discharge  Q, 
is  made  by  constructing  a line  having  the  direction  of  the 

solution  line  from  the  — curve  at  depth  d to  intersect 

a 

the  axis  of  ordinates.  The  value  read  at  the  intersection 
of  the  axis  of  ordinates  is  the  specific  energy  head  Hg^Q^. 


^Qd  = d + 


[v2] 

A i 

r q2  i 

_2g_ 

= 0.-1- 

_2ga2_ 

The  proof  of  this  construction  is  readily  seen  by  observing 
that  the  tangent  of  the  angle  formed  by  the  solution  line 


with  the  horizontal  is 


or 


2g 


the 


curve  at  d is 


1 


and  the  value  of  the  abscissa  of 


a 


Therefore, 


Q2" 

1 

[v2] 

[_2gJ 

* a2  ~ 

L2gJ 

Procedure  12. 


Generally,  when  the  specific  energy  head  H^q  and  the  dis- 
charge Q,  are  given,  two  possible  corresponding  depths  can  be 
ascertained.  The  graphical  determination  of  the  two  depths 
corresponding  to  a given  specific  energy  head  Hq^q  and  dis- 
charge Q is  given. 

a.  From  the  value  of  the  specific  energy  head  read  on 
the  axis  of  ordinates,  draw  a line  having  the  direc- 
tion of  the  solution  line  to  intersect  the  curve 

a2 

(usually  two  points). 

b.  Read  the  corresponding  depths  at  the  intersections  of 

the  — curve, 
a 


0 
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Procedure  13. 


The  corresponding  discharge  Q can  be  graphically  determined 
when  the  specific  energy  head  He,d  ami  depth  .d  are  given. 

Determine  the  discharge  Q corresponding  to  the  direction  of 
the  line  drawn  from  the  value  on  the  axis  of  ordinates 

to  the  point  on  the  ~ curve  at  depth  d.  (See  procedure  10 ). 

a 

Procedure  id. 


The  unique  situation  of  determining  the  critical  discharge  and 
critical  depth  corresponding  to  a given  specific  energy  head 
frequently  arises  in  hydraulic  computations.  An  example  of 
the  need  for  such  a computation  is  the  determination  of  the  dis- 
charge for  a reservoir  through  a control  section  into  a channel. 
A direct  solution  for  the  critical  discharge  and  depth  is  made 
graphically  by  the  use  of  ES-charts. 


The  graphical  solution  for  the  critical  discharge  Qc  and 
depth  dc  corresponding  to  the  specific  energy  head  He  is 
made  by  mentally  redesignating  the  axis  of  ordinates  as 
values  of  specific  energy  head  and  following  the  proce- 
dure below. 


a.  From  the  value  of  the  specific  energy  head  on  the 

axis  of  ordinates,  draw  a line  tangent  to  the 
curve . a 

b.  Determine  the  discharge  Qc  corresponding  to  the 
direction  of  the  solution  line  parallel  to  the 
tangent.  (See  procedure  10.) 

c.  Read  the  value  of  dc  corresponding  to  the  dis- 
charge Qc  determined  in  b.  (See  procedure  2). 

Though  it  is  noted  that  the  solution  line  drawn  in  step  a 

is  tangent  to  the  ~ curve  at  depth  dc,  the  point  of  tan- 
a 

gency  is  not  well  defined.  Solution  by  the  given  procedure 
is  the  most  accurate. 


Procedure  IQ. 


The  graphical  solution  for  the  velocity  head 


of  a given 


channel  corresponding  to  the  depth  d and  discharge  Q is  made 
by  determining  the  specific  energy  head  (see  procedure  11 ) 
and  subtracting  the  depth  d. 


Procedure  16. 


The  determination  of  a critical  slope  corresponding  to  a given 
Manning’s  coefficient  n and  discharge  Q (or  depth  d)  is  always 
unique . The  determination  of  the  discharge  Q,  (or  depth  d)  for 
which  a given  sQ  is  the  critical  slope  may  have  no  solution,  one 


solution,  or  two  solutions . 


The  determination  of  the  discharges  Q,  (or  depths  d)  for 
which  the  given  bottom  slope  sD  is  critical  is,  arithmeti- 
cally, a trial-and-error  solution.  This  procedure  gives 


a graphical  solution  for  any  possible  solutions  for 
either  discharge  Q,  or  depth  d within  the  depth  range 
the  charts. 

of 

a. 

Arbitrarily  select  values  of  depth  d for  the 
range  of  the  chart.  Prepare  a tabulation  of 
^c,  &>  Qn,  d anc^  [ Qn,  d ~ Qc,  dJ  - 

depth 

a. 

b. 

Plot  values  of  ([Qn,  d ~ Qc,dl  vs  d.)  . 

c . 

Read  the  solutions  of  depths  d for  which  the 
bottom  slope  sQ  is  critical  at  [Qn^q  - QCjq] 

given 
= 0 

from  the  chart  prepared  in  step  b.  The  discharges  Q 
for  which  sQ  is  a critical  slope  are  Qc?q. 

Water  Surface  Profile  in  a Computational  Reach 

Water  surface  profiles  for  a given  discharge  Q are  frequently  required 
in  the  work  of  the  Soil  Conservation  Service . They  are  required  in 
drainage,  irrigation,  flood  prevention,  and  other  conservation  works 
to  determine  depths  of  flow.  Water  surface  profiles  are  used: 

1.  to  ascertain  the  capacity  of  a channel, 

2.  to  ascertain  tailwater  depth  on  a hydraulic  structure, 

3-  to  ascertain  actual  tractive  forces  on  the  bottom  and  sides 
of  a channel, 

k.  to  ascertain  the  discharge  in  a delivery  channel, 

5.  to  locate  a hydraulic  jump,  and 

6.  for  many  other  objectives. 

The  water  surface  profile  is  determined  for  a steady  discharge  Q. 

Thus,  for  every  steady  discharge  there  is  a corresponding  water  surface 
profile.  Further,  a water  surface  profile  depends  on  a starting  depth. 
A starting  depth  is  the  depth  of  flow  which  is  given  or  pre -determined 
before  profile  computations  are  started.  It  is  emphasized  again  that 
the  water  surface  profile  depends  on  the  discharge  and  the  starting 
depth . 

The  determination  of  a water  surface  profile  is  merely  the  determina- 
tion of  the  depth  of  flow  at  selected  channel  sections  when  the 
channel  is  conveying  the  discharge  Q.  Two  consecutive  selected  sec- 
tions form  a computational  reach.  When  the  starting  depth  of  flow  is 
known  (or  has  been  determined)  at  one  of  these  two  sections  the  depth 
is  determined  at  the  other  section.  This  determined  depth  then  be- 
comes a known  depth,  which  is  used  as  a starting  depth  for  the  next 
computational  reach,  and  the  depth  is  determined  at  the  next  selected 
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section.  This  is  repeated  for  each  computational  reach  until  the 
depth  has  been  determined  at  every  selected  section. 

The  method  of  determining  the  water  surface  profile  in  a channel  when 
the  starting  depth  is  not  known  is  given  by  Procedure  l8,  which  re- 
quires a complete  understanding  of  Procedure  17. 

Procedure  17  is  used  to  determine  the  unknown  depth  at  the  end  of  the 
computational  reach  when  the  starting  depth  is  known.  The  determina- 
tion of  the  unknown  depth  involves  an  integration.  Since  it  is  not 
possible  to  formally  integrate  the  general  differential  equation  of 
varied  flow,  the  solution  is  often  obtained  by  numerical  integration. 

The  numerical  integration  is  an  approximation  that  may  be  made  as  accu- 
rate as  desired  by  selecting  computational  reaches  "sufficiently  short". 
The  charts  ES-125  through  ES-134  are  designed  to  give  a graphical  solu- 
tion for  this  integration.  The  graphical  solution  may  be  made  as  accu- 
rate as  the  charts  can  be  read,  assuming  the  computational  reach  {i2-lx) 
is  "sufficiently  short".  Whether  or  not  a reach  is  "sufficiently  short" 
depends  on  the  degree  of  accuracy  desired  by  the  user.  The  process  of 
determining  the  depth  by  the  use  of  the  chart  also  furnishes  the  friction 
losses  to  indicate  the  degree  of  accuracy  of  the  depth  determination. 

This  will  be  considered  after  the  explanation  of  the  depth  determination. 

The  graphical  determination  of  a depth  of  flow  at  the  end  of  a "suffi- 
ciently short  reach"  requires  the  following  data  to  be  given: 

1.  the  steady  discharge  Q - cfs, 

2.  a starting  depth  d - ft, 

3.  cross-sectional  shape  and  size, 

4.  Manning's  coefficient  n, 

5.  bottom  slope  of  the  channel  sQ  - ft/ft,  and 

6.  the  tolerable  error  to  be  permitted  in  the  depth  determination. 
Procedure  17. 

The  procedure  for  determining  the  depth  when  the  starting  depth 
is  known  is  given  on  sheet  3 of  ES-142.  It  assumes  the  user 
selects  the  value  of  n2(i2-i1)  which  in  turn  fixes  the  value  of 
the  length  of  the  computational  reach  {lz-lx)  i.e.,  the  distance 
between  the  stations  containing  the  starting  depth  and  the  depth 
to  be  determined. 

When  the  starting  depth  is  less  than  dc^Q  flow  is  super-critical 
and  computations  must  be  carried  in  a downstream  direction.  In 
other  words,  the  section  containing  the  depth  to  be  determined 
is  downstream  a distance  of  from  the  section  containing 

the  starting  depth.  When  the  starting  depth  is  greater  than 
dc^Q  flow  is  subcritical  and  computations  must  be  carried  in  an 
upstream  direction. 

drawing  the  solution  line  to  determine  the  depth  at  the  end  of 
the  computational  reach  also  furnishes  the  change  in  velocity 
head  and  the  friction  loss  in  the  reach.  See  sheet  1 of 
ES-142  which  illustrates  the  determination  of  velocity  heads 
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and.  friction  losses  when  flow  is  subcritical,  supernormal, 
accelerated,  and  on  a positive  sloped  channel.  The  graphical 
integration  of  energy  loss  "by  friction  is  an  approximation. 

It  depends  on  the  depth  of  flow  which  generally  varies 
throughout  the  reach  of  the  channel.  The  rate  of  energy  loss 
between  the  sections  of  the  starting  depth  and  the  section  of 
the  determined  depth  is  taken  as  the  arithmetical  average  of 
the  instantaneous  rate  of  friction  loss  of  the  two  end  sec- 
tions of  the  computational  reach  (see  procedure  8).  This  is 
equivalent  to  saying  that  the  friction  loss  for  half  of  the 
computational  reach  is  evaluated  with  a depth  of  flow  equal 
to  the  starting  depth,  and  the  friction  loss  for  the  other 
half  of  the  computational  reach  is  evaluated  with  a depth  of 
flow  equal  to  the  determined  depth  of  flow. 


The  error  in  this  approximation  is  always  less  than  the  differ- 
ence of  these  two  evaluated  friction  losses  Ih-p  - hf  I,  since 

1 2 1 1 


the  depth  varies  monotonically . 


The  instantaneous  rate  of  friction  loss  for  the  depth  of  flow 
d is 


dhp 

dJ~ 


s0 


Q 

Qn,  d_ 


2 


The  friction  loss  for  the  distance 


hf 


r q i 

2 

So 

_Qn,d_ 

2 

The  distance  between  the 


curve) is 


_<£,d 


g^o-'O 


at 


-curve  and  the  Up-curve  (or  Up- 
the  depth  d,  and  the  tangent  of 


the  angle  formed  by  the  solution  line  with  the  horizontal  is 
. The  friction  loss  for  a depth  d over  the  distance 


l_2g. 


’(VO 


is  the  vertical  component  of  the  projection  on  the 


solution  line  of  the  horizontal  line  segment  between  the  — ~ - 

a 


curve  and  the  Up-curve  (or  the  — -curve  and  the  Up-curve). 

ci 

The  energy  relations  of  varied  flow  shown  on  sheet  1 of 
ES-1^-2  is  for  one  case  of  the  twelve  types  of  water  surface 


5-6 


curves  in  prismatic  channels.  Sheet  2 of  ES-1^2  summarizes  the  typical 
graphical  procedures  for  depth  determination  for  all  cases  of  varied 
flow.  Observe  that  for  subcritical  flow  the  starting  point  A is  on  the 
Up-curves,  and  for  supercritical  flow  the  starting  point  A is  on  the 
Up-curves . 


Water  Surface  Profiles  in  a System  of  Channels 

The  water  surface  profile  is  often  desired  for  a system  of  channels 
which  may  contain  hydraulic  structures  and  pools.  The  solution  for 
profiles  for  such  a system  of  channels  resolves  into  one  of  three  prob- 
lems, that  is: 

1.  The  determination  of  the  water  surface  profile  corresponding 
to  a given  steady  discharge  Q,  and  a tailwater  depth  when  flow 
is  subcritical,  or  the  determination  of  a water  surface  pro- 
file corresponding  to  a given  discharge  Q and  headwater  depth 
when  flow  is  supercritical. 

2 . The  determination  of  the  steady  discharge  Q,  in  a channel  for 
a headwater  depth  and  a tailwater  depth  corresponding  to  the 
unknown  discharge  Q.  This  is  often  a trial-and-error  problem 
which  involves  only  the  principles  of  the  first  stated  prob- 
lem for  selected  discharges.  Examples  of  problems  which  have 
a given  headwater  and  tailwater  depth  instead  of  a given  dis- 
charge are : 

a.  When  flow  is  subcritical  in  a bifurcated  channel  a 
headwater  condition  is  imposed  at  the  section  of  fork- 
ing. The  headwater  condition  at  the  bifurcated  section 
is  usually  assumed  to  be  level.  This  fact  is  used  to 
confirm  the  headwater  depths  and  determine  the  divided 
discharges,  which  also  depend  on  the  tailwater  depth. 
This  is  a water  surface  profile  problem. 

b . A delivery  channel  from  a pool  has  a given  headwater 
elevation  at  the  outlet  of  the  pool.  The  elevation  of 
the  water  surface  in  the  pool  along  with  the  tailwater 
depth  will  determine  the  discharge  in  the  channel. 

When  the  outlet  of  the  pool  is  not  a control  section, 
this  becomes  a water  surface  profile  problem.  It  was 
tacitly  assumed  no  outlet  structure  exists  at  the  en- 
trance of  the  channel  in  this  example. 

The  determination  of  a tailwater  depth  for  a given  discharge  Q 
and  headwater  depth  when  flow  is  subcritical,  or  the  determina- 
tion of  a headwater  depth  for  a given  discharge  Q and  tailwater 
depth  when  flow  is  supercritical.  For  prismatic  channels  hav- 
ing no  break  in  grade,  this  solution  can  be  obtained  by  one 
water  surface  profile  determination. 
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Procedure  l8. 


For  profile  determination,  a channel  containing  hydraulic 
structures  and/or  pools  is  divided  into  reaches  which  do  not 
contain  a structure  nor  a pool.  These  reaches  are  called 
channels  under  consideration  and  are  considered  individually 
and  successively  in  an  upstream  direction.  A channel  under 
consideration  may  or  may  not  have  a structure  or  pool  at  its 
end  section.  The  water  surface  profile  depends  on  the  dis- 
charge and  starting  depth.  Since  structures  and  pools  at  the 
ends  of  a channel  under  consideration  can  affect  the  discharge 
and  starting  depth,  these  effects  are  evaluated  before  commenc- 
ing water  surface  profile  computations.  A brief  description 
of  the  effects  of  pools  and  the  various  types  of  structures 
on  discharge  and  starting  depths  are  considered. 

The  evaluation  of  the  pool  level  will  determine  the  discharge 
Q in  the  channel  conveying  or  delivering  water  from  the  pool. 

The  determination  of  the  discharge  Q from  the  pool  is  simple 
when  the  entrance  to  the  channel  is  known  to  be  a control 
section.  For  this  situation  the  pool  depth  at  the  entrance 
of  the  delivery  is  H6^q.  When  the  entrance  section  is  not 
a control  section,  a rating  curve  is  computed  for  the  entrance 
section.  This  requires  that  water  surface  profiles  correspond- 
ing to  various  discharges  be  computed  upstream  to  the  pool. 

The  overfall  type  of  hydraulic  structure  in  a channel,  whether 
located  at  the  outlet  of  a pool  or  not,  often  is  a control 
section  at  which  the  depth  of  flow  at  the  weir  is  known  to  be 
critical  depth.  If  the  elevation  of  the  tailwater  on  the 
overfall  type  of  structure  is  greater  than  the  elevation  of 
the  critical  depth  over  the  weir,  the  weir  section  is  not  a 
control  section;  and  the  water  surface  profile  corresponding 
to  various  discharges  are  computed  upstream  to  the  weir  section. 
A rating  curve  for  the  submerged  weir  section  is  computed 
from  these  profiles. 

Sluice  structures  in  a channel  usually  control  the  discharge 
Q and  the  starting  depth  at  the  upstream  end  of  a channel 
under  consideration.  When  they  do  not  control  the  starting 
depth,  the  tailwater  on  the  sluice  must  be  computed  to  deter- 
mine the  discharge.  Again,  this  is  accomplished  by  computing 
water  surface  profiles  to  obtain  the  tailwater  on  the  sluice. 

The  headwater  on  the  sluice  can  then  be  determined  by  its 
hydraulic  characteristics  for  the  given  discharge. 

Closed  conduits  may  or  may  not  flow  full.  Water  surface  pro- 
files may  need  to  be  computed  in  partially  filled  closed  con- 
duits to  determine  rating  curves.  Technical.  Release  No.  14 
gives  the  method  of  obtaining  the  rating  curve  for  rectangular 
and  circular  culverts  with  all  possible  tailwater  conditions 
by  the  use  of  a digital  computer.  When  the  closed  conduit  is 
flowing  full  and  the  tailwater  has  no  control  on  the  discharge, 
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the  discharge  is  determined  entirely  by  the  hydraulic  charac- 
teristics of  the  structure.  When  the  tailwater  does  have  an 
effect  on  the  discharge,  then  the  water  surface  profile  com- 
putations are  required  to  obtain  the  discharge  of  the  struc- 
ture . 

When  pools  and  hydraulic  structures  solely  control  the  dis- 
charge, the  profile  for  the  channels  under  consideration 
connecting  the  pool  or  structure  is  determined  independently. 
However,  it  is  not  always  possible  to  ascertain  before  the 
profile  computations  are  made  whether  or  not  the  pool  or  struc- 
ture solely  controls  the  discharge. 

The  same  given  data  is  required  for  water  surface  profile  compu- 
tations in  each  of  the  channels  under  consideration  as  is  re- 
quired for  computational  reaches. 

This  data  is: 

1.  the  steady  discharge  Q - cfs, 

2.  a starting  depth  d - ft, 

J.  cross-sectional  shape  and  size, 

4.  Manning’s  coefficient  n, 

5.  bottom  slope  of  the  channel,  and 

6.  the  tolerable  error  to  be  permitted  in  the  depth  deter- 
mination. 

Often  the  discharge  Q is  not  given  or  the  starting  depth  is  not 
given.  The  following  are  the  three  possible  cases. 

a.  When  the  discharge  Q is  given  and  the  starting 
depth  is  not  given  (that  is  a tailwater  or  head- 
water depth)  and  no  structures  are  involved  at 
the  starting  end  of  the  channel  under  considera- 
tion, the  starting  depth  may  be  determined  by 
the  procedure  given  on  sheets  4,  5 > and  6 of 
ES-142.  In  this  situation  it  is  necessary  to 
know  the  tailwater  condition  when  flow  is  sub- 
critical  and  the  headwater  condition  when  flow 

is  supercritical  to  determine  the  starting  depth. 
This  requires  that  sufficient  length  of  channel 
beyond  the  channel  under  consideration  be  given. 
When  the  downstream  end  section  of  the  channel 
under  consideration  having  subcritical  flow  ter- 
minates at  a structure  or  pool,  the  starting 
depth  can  be  determined  by  the  hydraulic  charac- 
teristics of  the  structure  or  the  elevation  of 
the  pool. 

b.  When  the  discharge  Q is  not  given  and  the  start- 
ing depth  is  given,  the  discharge  Q is  determin- 
able in  subcritical  flow  if  the  headwater  depth 
(or  some  other  elevation  condition)  is  given. 
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If  the  starting  depth  is  given  in  supercritical  flow, 
the  tailwater  depth  (or  some  other  elevation  condition) 
must  also  he  given  to  determine  the  discharge  Q.  Select 
various  discharges,  and  compute  their  corresponding  water 
surface  profiles  from  the  given  starting  depth.  A rating 
curve  is  prepared  at  the  station  of  the  given  elevation 
condition  to  determine  the  discharge  Q.  Sometimes  the  dis- 
charge is  controlled  by  a structure  or  by  a control  section 
at  the  upstream  end  of  the  channel  under  consideration  for 
the  given  "elevation  condition",  and  the  tailwater  has  no 
effect  on  the  discharge.  In  these  situations  the  discharge 
is  determined  without  recourse  to  water  surface  profile 
determination. 

c.  When  the  discharge  Q,  and  the  starting  depth  are  not  given, 
an  elevation  condition  other  than  what  would  normally  be 
used  for  a starting  depth  must  be  given.  The  discharge  Q 
and  starting  depth  can  be  determined  for  only  those  cases 
in  which  no  hydraulic  jump  occurs  between  the  starting 
depth  to  be  determined  and  the  known  elevation  condition. 
Various  discharges  are  selected,  and  their  starting  depths 
are  determined  by  the  procedure  given  on  sheets  4,  5 y and 
6 of  ES-142.  The  corresponding  water  surface  profiles  for 
the  selected  discharges  and  their  corresponding  starting 
depths  are  computed.  A rating  curve  is  prepared  at  the 
station  of  the  given  elevation  condition  to  determine  the 
discharge  Q,. 


« 


3-u 


HYDRAULICS1  n2(l2-A)  VS.  (j22-W  for  various  values  of  n 
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HYDRAULICS:  Energy  relations  of  varied  flow 
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HYDRAULICS : Summary  of  graphic  procedures  on  Water  Surface  Profile  charts  ( ES- 125  thru  ES-134)  for  prismatic  channels 

Subcritical  Flow  Supercritical  Flow 
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HYDRAULICS:  Procedure  for  computing  o water  surface  profile  for  a channel  under  consideration  when  the  starting  depth  is  or  is  not  given. 
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chap™  4 - EXAMPLES 


EXAMPLE  1 

Given:  (l)  Trapezoidal  channel,  z = 2 

(2)  b = 50  ft 
(5)  d = 5-2  ft 

Determine : The  critical  discharge  Qc ^ ^ corresponding  to  the  given 

depth . 

Solution:  Select  the  proper  chart  (ES-129,  sheet  140).  On  the  left- 

hand  line  chart  (Qc  £ vs  d)  for  a value  of  d = 5*2  ft,  read 
Qc, d = !7^0  cfs. 

EXAMP1PI  2 

Given:  (l)  Trapezoidal  channel,  z = 2 

(2)  b = 50  ft 
(5)  Q = 1700  cfs 

Determine:  The  critical  depth  dc  q corresponding  to  the  given  dis- 

charge . 

Solution:  Select  the  proper  chart  (ES-129,  sheet  140).  Mentally 

redesignate  the  parameters  (Qc  ^ vs  d)  of  the  left-hand  line  chart 
as  (Q  vs  d^q).  For  Q = 1700  cfs,  read  dc^Q  = 5-16  ft. 

EXAMPLE  5 

Given:  (l)  Trapezoidal  channel,  z = 2 

(2)  b = 50  ft 

(5)  d = 5.9  ft 

(4)  sQ  = 0.005  ft/ft 

(5)  n = 0.055 

Detemine:  The  normal  discharge  Qn  ^ corresponding  to  the  given 

depth . 


Solution:  Select  the  proper  chart  (ES-129,  sheet  140).  On  the  right- 


hand  line  chart  d vs 


( 

d 

\ 

nQn,  d 

d vs 

t±l 

j for  d = 5*9  ft,  read 

V 

w 1 

H 

0 

Lso1/2 * 4 5] 

= 750. 


Obtain 


n 


s01/2 


=0.64  from  ES-158. 

n^n,d 


s 1/2 
bO 


n 


so 1/2 


750 

0.64 


= 1172  cfs. 


Solve  for  Qn^  ^ = 
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EXAMPLE  4 

Given:  (l)  Trapezoidal  channel,  z = 2 

(2)  b = 50  ft 

(3)  s0  = 0.003  ft/ ft 

(4)  Q = 1500  cfs 

(5)  n = 0.035 

Determine : The  normal  depth  dn?Q  corresponding  to  the  given  discharge, 


Solution:  From  ES-I38  find 


Compute 


n 


s x/2 
Lso 


- 0.64 


nQ 

n 

so1/2. 

1 — 1 

0 

Q = (0.64) (1500)  = 960. 


Select  the  proper  chart  (ES-129,  sheet  140). 

~p  ^,d 


Mentally  redesignate  the  parameters  ^d  vs 

nQ 


Sr,1/2 


L °o 


hand 

dn,Q 

line 

= 4. 

chart  as  | dn  q vs  - 

v ' L 

50  ft. 

EXAMPLE  5 

Given : 

(1) 

Trapezoidal  channel. 

(2) 

b = 50  ft 

(5) 

n = 0.04 

00 

d = 4.16  ft 

(5) 

Q = 1700  cfs 

For 


nQ 


s01/2 


of  the  right- 
= 960,  read 


Determine : The  normal  slope  sn  corresponding  to  the  given  discharge 

and  depth. 


Solution:  Locate  the  proper  chart  (ES-129*  sheet  140). 

~n  Qn,d 


Mentally  redesignate  the  parameters  [ d vs 

nQ 


of  the  right- 


hand  line  chart  as  d vs 


.( sn,Qd)1/2_ 


s V2 

s0  ' 

For  d = 4.16  ft,  read 


nQ 


= 835. 


(sn,Qd)1/2_ 

(Method  2.)  Compute 


n 

nQ 

_( sn,  Qd)  1/2_ 

_(sn,  Qd)1/2_ 

1 _ _ 


Q 1700 


= 0.491 
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3 


EL 


Mentally  redesignate  the  parameters  [ sQ,  n,  and 

as  (sn,Qd>  n>  -4- P'  JV  For 

\ U sn,  Qd'  jJ 

n = 0.04,  read  = 0.0066  ft/ft. 


n 


So1/2 


tof  ES-I38 


>n,Qd)1/2 


3 


= 0.491  and 


EXAMPLE  6 

Given:  (l)  Trapezoidal  channel,  z = 2 

(2)  b = 50  ft 

(3)  n = 0.017 

(4)  d = 4.53  ft 

Determine : The  critical  slope  sc^ ^ corresponding  to  the  given  depth. 

Solution:  Select  the  proper  chart.  (ES-129,  sheet  140). 

On  the  left-hand  line  chart  (Qc  ^ vs  d)  for  d = 4.53  ft,  read 
Qc,d  = 3000  cfs. 

n Qn,  d 


Mentally  redesignate  the  parameters  f d vs 


Sr,1/2 


jof  the  right-hand 


line  chart  as  d vs 


n QC,d 


L(sc,d)1/2 


n Q, 


c,  d 


Led) 1/2 


= 970, 


For  d = 4.53  ft,  read 

( Method  2 . ) Compute 

Mentally  redesignate  the  parameters  ^s0,  n>  and. 
as  I sc  n,  and  

V ’ L<sc, 


n 

n Qc,d‘ 

1 970 

(Bed)1/2. 

(Bed)1/2. 

Qc, d 3000 

n 


s1/2 


= 0.323 


of  ES-138 


n 


d)1/2 


For 


n 


(Bed)1/2 


= 0.325  and  n = 0.017,  read  sc?a  = 0.00277  ft/ft. 


J 


EXAMPLE  7 


Given:  (l) 

Trapezoidal  channel 

(2) 

b = 50  ft 

(3) 

n = 0.017 

• (4) 

Q = 3000  cfs 

4-1+ 


Determine : The  critical  slope  sCjq  corresponding  to  the  given 

discharge. 

Solution:  Select  the  proper  chart.  (ES-129,  sheet  l40). 

Mentally  redesignate  the  parameters  (Qc  ^ vs  d)  of  the  left-hand 
line  chart  as  (Q  vs  dCjQ).  For  Q = 3000  cfs,  read  dc^Q  = 4.53  ft. 

n Qn,d~  ’ 


Mentally  redesignate  the  parameters  [ d vs 

nQ 


s x/z 

L sO  _ 


line  chart  as  jdc?Q  vs 

\ 

nQ. 


of  the  right-hand 
For  dc^Q  = 4.55  ft,  read 


= 970. 


(sc,  q)1/2_ 

(Method  2.)  Compute 
Mentally  redesignate  the  parameters  ( s0,  n,  and 


n 

nQ 

1 970 

_(sc  ,Q)1/2_ 

Jscq)1/2^ 

Q 3000 

n 


s V2 
LrO 


lof  ES-I38  as 


fsc.Q,  n,  and  — — — — — ).  For  

\ N,q)172J/  |_(sc,q) 

read  sCjq  = 0.00277  ft/ft. 

EXAMPLE  8 

Given : (l)  Trapezoidal  channel,  z = 2 

(2)  b = 50  ft 

(3)  n = 0.020 

(4)  d = 4.50  ft 

(5)  Q = 2000  cfs 

dhf 


1/2 


= 0.323  and  n = 0.017, 


Determine : The  rate  of  friction  loss 

depth  and  discharge. 


di 


corresponding  to  the  given 


Solution:  Locate  the  proper  sheet  (ES-129,  sheet  l40).  On  the  right- 


hand  line  chart 


n(4n,d 


vs  d for  d = 4.5  ft,  read 


nQn,  d 
so 


1/2 


= 957- 


Qn,d 


Compute 


Compute  s = sc 


nQn,d 
■so 1/2 


997 

0.020 


= 2,290,000,000. 


_Qn,d. 


or 


<d 


QQQ;  QQ.Q 0 001747  ft /ft 

2,290,000,000  ~ °-oolw 
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EXAMPLE  9 

Given:  (l) 

Trapezoidal  channel 

(2) 

b = 50  ft 

(3) 

Q = 2100  cfs 

Determine : The  direction  of  the  solution  line  for  the  given  discharge. 

Solution:  Select  the  proper  sheet  (ES-129>  sheet  140).  The  smallest 

value  of  m is  5.0  Compute  — = — 7^  = 700  which  can  he  read  on  the 

0 m ^ 

^ scale.  The  direction  of  the  solution  line  corresponding  to  the 
discharge  2100  cfs  is  parallel  to  the  vector  drawn  from  the  point 
m = 5 to  the  point  ^ = 700. 


EXAMPLE  10 

Given : (l)  Trapezoidal  channel,  z = 2 

(2)  h = 50  ft 

(5)  Direction  of  the  solution  line  defined  hy  the  two  points: 

a.  ■—  curve  at  d = 4 ft  and 
ar 

h.  axis  of  ordinates  at  the  value  of  7 ft. 


Determine : The  discharge  Q corresponding  to  the  given  direction  of 

the  solution  line. 


Solution:  Select  the  proper  chart  (ES-129>  sheet  140). 

Locate  the  two  given  points  which  define  the  given  direction  of  the 
solution  line. 

, Q 

Draw  a line  showing  this  direction  from  m = 4 to  intersect  the  — 

G 0 m 

scale  at  the  value  ^ = 805«  Multiply  the  value  — = 805  hy 

m = 4.0  to  obtain  the  corresponding  discharge  Q = 5220  cfs. 


EXAMPLE  11 


Given:  (l) 

Trapezoidal  channel 

(2) 

b = 50  ft 

(3) 

Q = 2100  cfs 

(4) 

d = 5.8  ft 

Determine : The  specific  energy  head  He  corresponding  to  the 

given  depth  and  discharge. 

Solution:  Select  the  proper  chart  (ES-129^  sheet  l40).  Determine 

the  direction  of  the  solution  line  for  the  discharge  Q = 2100  cfs. 
(See  example  9* ) Construct  a line  having  the  direction  of  the 

solution  line  from  the  ■—  curve  at  d = 5*8  to  intersect  the  axis 

ar 

of  ordinates.  The  value  read  at  this  point  of  intersection  is 
the  specific  energy  head,  He  = 5-24  feet. 
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EXAMPLE  12 


Given:  (l) 

Trapezoidal  channel 

(2) 

b = 50  ft 

(3) 

He,Q  = 5-6  ft 

(4) 

Q = 2100  cfs 

Determine : The  depths  d at  which  the  given  discharge  may  occur  in 

the  given  channel  at  the  given  specific  energy  head  Hs^q. 

Solution : Select  the  proper  chart  (ES-129;  sheet  l40) . 

Determine  the  direction  of  the  solution  line  for  the  given  dis- 
charge . ( See  example  9 * ) Draw  a line  having  the  direction  of 
the  solution  line  from  the  value  H6jq  = ^>.6  ft  on  the  axis  of  ordi- 
nates to  intersect  the  curve  at  two  points.  The  depths  read  at 
the  points  of  intersection  are: 

d = 2.85  ft  and 
d = 4.74  ft. 

EXAiffLE_l^ 

Given:  (l)  Trapezoidal  channel,  z = 2 

(2)  b = 50  ft 

(3)  d - 4.74  ft 

(4)  He,d  = 5.6  ft 

Determine : The  discharge  Q corresponding  to  the  given  depth  and 

the  given  specific  energy  head. 

Solution:  Select  the  proper  chart  (ES-129*  sheet  140). 

Draw  the  line  from  the  value  of  5*6  ft  on  the  axis  of  ordinates 

to  the  curve  at  the  depth  4.74  ft.  Determine  the  discharge 
a  *  2 3 

Q corresponding  to  the  direction  of  this  line,  (See  procedure  10.) 

Q 

On  selecting  in  = the  value  of  — is  equal  to  700.  Hence,  the 
corresponding  discharge  Q = J00  x ) = 2100  cfs. 

EXAMPLE  14 

Given:  (l)  Trapezoidal  channel,  z = 2 

(2)  b = 50  ft 

(3)  He  = 5-5  ft 

Determine : The  critical  discharge  Qc  and  the  critical  depth  dc 

corresponding  to  the  specific  energy  head  Hg  = 5*50  ft. 

Solution:  Select  the  proper  chart  (ES-129*  sheet  140).  Mentally 

redesignate  the  axis  of  ordinates  of  the  central  chart  d as  Hg. 

From  the  point  He  = 5*50  ft  on  the  axis  of  ordinates,  draw  the 

curve.  Determine  Qc  for  this  given 
a 


line  tangent  to  the 
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direction  of  the  solution  line.  Qc  = 2286  cfs.  (See  procedure  10o) 

Mentally  redesignate  the  parameters  (d  vs.  Qc  d)  of  the  left-hand 
line  chart  as  (dc  vs.  Qc ) . For  Qc  = 2286  cfs,  read  dc  = 3-82  ft. 
(See  procedure  2.) 


) 


) 


EXAMPLE  13 
Given: 


(1)  Trapezoidal  channel,  z = 2 

(2)  b = 50  ft 

(3)  Q = 2100  cfs 

(4)  d = 3.8  ft 


Determine : The  velocity  head 


fv2‘ 


_2&. 


of  the  given  channel  corresponding 


to  the  given  depth  and  discharge. 


Solution:  Select  the  proper  chart  (ES-129,  sheet  140). 

Determine  the  specific  energy  head  Hg  ^ corresponding  to  the  given 

depth  and  discharge.  (See  example  11.) 


He,Qd  = d + 


5.24  ft 


"v£ 

2g_ 


He,Qd  - d = 5-24  - 3-80  = 1.44  ft 


EXAMPLE  16 


Given:  (l) 

Trapezoidal  channel,  z = 2 

(2) 

b = 50  ft 

(5) 

n =0.020 

(4) 

sQ  = 0.0025  ft/ft 

Determine : The  depth  d and  discharge  Q corresponding  to  the  critical 

slope  equal  to  sQ . 


Solution:  Select  the  proper  charts  (ES-129*  sheets  139  through  144). 

Determine  Qc,d  arL(3-  Qn,  d for  several  values  of  d,  using  a tabular 

form.  Column  2 is  read  directly  from  the  left-hand  line  chart, 
and  Column  3 is  read  directly  from  the  right-hand  line  chart. 


Compute 


n 

r 0.020  1 

q x/2 

Lso  ' _j 

_(0.0025)1/2_ 

= 0.40 


Column  4 is  Column  3 divided  by  the  computed  value  of 


n 


s 1/2 
so 


Column  5 is  the  difference  of  Qc^d  and  Qn^d  for  each  depth  d. 


Since  these  differences  are  small  in  relation  to  the  total  dis- 
charges, the  diff  ^rences  are  plotted  against  depth  to  find  a more 
accurate  intersection  of  the  points  where  dc^Q  = dn^Q. 
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( 


(1)  (2)  (3)  W (5) 


a 


Qc,a 


"n  Qn,a~ 

s 1/2 
so 


On,  a ( On,  a - Qc,a  ) 


29 

77,200 

31,  300 

78,250 

+ 1050 

25 

57,900 

23, 100 

57, 750 

- 150 

20 

38, 000 

17,  750 

36,875 

- 1125 

15 

22, 500 

8,720 

21, 050 

- 1750 

10 

11, 100 

3,950 

9,875 

- 1225 

5 

3, 525 

1,150 

2,875 

- 650 

3 

1,  570 

776 

1, 190 

- 380 

2 

837 

272 

605 

- 252 

1 

289 

75 

187 

- 102 

•5 

101 

23-5 

59 

72 

.2 

25.8 

5 

12.5 

13 

Plot  the  curve  ((Qn, a “ Qc,d)  vs.  dj  . 


Read  the  depth  a = 25.6  ft  corresponaing  to  the  critical  slope  equal 
to  sQ  from  the  plottea  curve  at  (Qn^ a ~ Qc,d)  = 0.  Read  the  discharge 
Q corresponding  to  the  depth  of  25.6  ft  on  the  given  slope  directly 
from  the  left-hand  line  chart. 

The  discharge  Q,  for  which  the  given  bottom  slope  sQ  is  critical,  is 
Q = 60, 600  cfs . 


( 
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1 


) 


EXAMPLE  17 


Given: 


(1) 

(2) 

(3 

00 

(5) 

(6) 
(7) 


Trapezoidal  channel,  z = 2 

b = 50  ft 

sQ  = 0,005  ft/ft 

Q = 1825  cfs 

n = 0.040 

d = 4.15  ft 

Tolerable  error  of  0.05  ft 


Determine : The  unknown  depth  at  the  distance  from  the 

section  of  known  depth. 


Solution:  Select  the  proper  chart  (ES-129*  sheet  140). 

Select  n2(i2-f1)  = 0.10  From  ES-I37,  (i  -i1)  = 62.5  ft. 

Compute  s0(i2-i1)  = 0.005  (62.5)  = 0.512  ft. 

Read  dc^Q  corresponding  to  the  discharge  Q = 1825  cfs  on  the  left- 
hand  line  chart. 

dc,Q  = 3.30  ft  < d = 4.15  ft 

Thus,  flow  is  subcritical  and  water  surface  profile  computations 
are  made  in  an  upstream  direction. 


Find  point  A on  the  chart  at  d = 4.15  ft  on  the  Up-curve  for 
n2(i2-i1)  = 0.10  (See  page  4-10.) 

Scale  the  distance  s0(£2-£1)  = 0.^12  down  from  point  A.  This  is 
point  B,  at  d = 3*838  ft. 


Determine  the  direction  of  the  solution  line  through  points  m = 3 

on  the  m-scale  and  — = 608  on  the  — scale.  Draw  the  solution  line 

m m 

through  point  B to  intersect  the  conjugate  Up-curve 

(n 2(i2-i1)  = O.io)  . Read  the  depth  d-,_,  which  is  greater  than  dc^Q, 

as  4.39  ft.  This  is  the  depth  at  the  station  62.5  ft  upstream  from 
the  section  at  which  d2  is  4.15  ft. 

Project  the  horizontal  line  segment  between  the  -%■  curve  and  A to 

a 

the  solution  line.  The  vertical  component  of  this  projected  line 

segment  is  4.08  - 3*84  = 0.24  = h-f  . 

2 

Project  the  horizontal  line  segment  between  C and  the  -h?  curve  to 

cl 

the  solution  line.  The  vertical  component  of  this  projected  line 

segment  is  4.39  - 4.19  = 0.20  = h-p  . 

1 

This  solution  has  an  error  of  less  than  I h-p  - h-p  I = 

2 1 

| 0.24  - 0.20 | = 0.40  ft,  which  is  less  than  the  tolerable  error 
of  0.05  ft. 


12+50 
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EXAMPLE  18a 


Given: 


(1)  Trapezoidal  channel,  z = 2 

(2)  b =’45  ft 

(5)  The  profile  of  the  channel  bottom  as  shown  in  the  figure 

(4)  Q = 4000  cfs 

(5)  n = 0.020 

(6)  The  elevation  of  a lake,  into  which  the  channel  dis- 
charges, is  119.4  ft  at  approximately  six  miles  down- 
stream from  the  channel  under  consideration.  The 
channel  slopes  were  observed  to  be  generally  uniform 
all  the  way  to  the  lake. 

(7)  Tolerable  error  of  0.05  ft. 


Determine : The  water  surface  profile  for  the  channel  under  consider- 

ation. 


Solution:  Since  there  is  no  given  starting  depth,  the  procedure 

given  on  sheet  4 of  ES-142  is  used  in  the  solution  of  this  example. 

The  critical  slope  corresponding  to  the  discharge  Q = 4000  cfs  is 
determined  (see  procedure  6). 

sc,Q  = 0.00180  ft/ft  > s0  = 0.0007 

Thus,  there  is  a mild  slope  downstream  from  the  first  break  in 
grade  below  the  downstream  end  of  the  channel  under  consideration 
(the  reach  between  69  + 25  and  79  + 25). 

The  selected  charnel  reach  = ([69  + 25]  - [64  + 25])  = 501  ft 

contains  the  break  in  grade  at  station  69  + 25. 
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Compute  dn  q for  the  bottom  store  s ^ - n nnnv  ( , \ 

dn,Q  = 9.otVt.  Compute  the  wafer  pLii?!  ^T*1”  k) 

station  69  + 26  using  the  starting  depth  a T-  9 fr°m 

depth  d = 5.98  ft  at  station  64  + 25  (see  & 

m hr^Sni^t^d^i?-*04  for  ttati°“s  ^ + 26 

depths  give  the  elevation  at  these  iStlo^'^^t®^ 


a^h“Tat00^TioT:  lndiTs  that  the 

the  tarlwater  conditions  would  Influence  the  depth  of  now  ft 'stolon 

& + * s-r  Sprfn  ^ T'z/riP 

starting  depth  d = 6.02  ft  at  station  64  ^ 25.  glVeS  ^ 

- rya«.“*«s  t„tsh  sz?*»  *■ 


U9+50 


• if -13 

Compute  the  profile  to  the  upstream  end  of  the  channel  under  con- 
sideration. A tabulation  of  reach  lengths  (£2-£-,)  and  correspond- 
ing depths  are  given  below  for  each  step.  Since  there  are  no  steep 
slopes  upstream,  the  profile  for  the  channel  under  consideration  is 
complete. 


n1 2 * 4 5 6(£2-ii) 

(Jz-lx) 

2(vi) 

Depth 

Station 

0 

0 

0 

6.02 

64  + 25 

0.05 

125 

125 

6.90 

63  + 00 

0.10 

250 

375 

7.60 

60  + 50 

0.10 

250 

625 

8.02 

58  + 00 

0.10 

250 

875 

8.35 

55  + 50 

0.10 

25O 

1125 

8.58 

53  + 00 

0.10 

250 

1575 

8.80 

50  + 50 

0.10 

250 

1625 

9.00 

.48  + 00 

0.10 

250 

1875 

9.15 

45  + 50 

0.02 

50 

1925 

9.18 

45  + 00 

0.40 

1000 

2925 

9.58 

35  + 00 

0.40 

1000 

5925 

9.52 

25  + 00 

0.40 

1000 

4925 

9.65 

15  + 00 

0.10 

250 

5175 

9.64 
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EXAMPLE  l8b 


Given: 


(1)  Trapezoidal  channel,  z = 2 

(2)  b = 20.0  ft 

(2)  n = 0.04 

(4)  Tailwater  depth  at  station  40  + 00  = 8,50  ft 

(5)  Pool  elevation  at  upstream  end  of  charnel  under  consider 
at ion  is  105«94 

(6)  The  profile  of  the  channel  under  consideration  is  shown 
in  the  following  figure. 


4-l4 


Determine : The  steady  discharge  Q in  the  channel  under  consideration 

for  the  given  taiTwater  depth  d and  headwater  condition. 


Solution:  Select  several  discharges  such  that  the  computed  pool 

elevations  for  two  of  the  selected  discharges  will  bracket  the 
known  pool  elevation. 

Select,  in  turn,  Q = 1000  cfs,  Q = 1500  cfs,  and  Q = 2000  cfs.  Start 
all  water  surface  profile  computations  at  station  40  + 00  with 
d = 8.50  ft,  and  compute  the  profiles  for  the  selected  discharges  to 
station  25  + 00.  (See  sheet  4 of  ES-142). 


Compute  the  pool  elevation  by  adding  the  depth  d and  the  velocity 
1^21 

head 


_2gJ 


(corresponding  to  the  selected  discharges)  to  the  channel 
bottom  elevation  97*00  ft  at  station  25  + 00. 

Plot  the  following  rating  curve  (pool  elevation  vs.  discharge  Q). 


1000  1500  2000 
Channel  Discharge  Q,  cfs. 

Rating  curve  of  pool  at  exit  station  25  + 00 

Read  channel  discharge  Q,  = 1750  cfs  corresponding  to  the  given  pool 
elevation  105-9^ 


25+00 


EXAMPLE  l8c 

Given:  (l) 

Trapezoidal  channel,  z = 2 

(2) 

b = 20  ft 

(5) 

n = 0.04 

(M 

Pool  elevation  at  upstream  end  of  channel  under  con- 
sideration is  106.11 

(5) 

The  profile  of  the  channel  under  consideration  is 
shown  in  the  following  figure. 

(6) 

Tolerable  error  = 0.10  ft 

Pool  Elevation  106»11 


Determine : The  discharge  Q of  the  channel  and  the  channel  depth 
of  flow  at  station  25  + 00. 

Solution:  Select  the  discharge  Q = 1000  cfs.  Since  there  is  no 

starting  depth  given,  the  computational  procedure  will  begin  as 
given  on  sheet  4 of  ES-142.  Begin  the  water  surface  profile 
computation  at  station  40  + 00  for  Q = 1000  cfs  with  an  assumed 
depth  dn^Q  = 7*97  ft  (where  sq  = 0.001).  Compute  the  water  sur- 
face profile  upstream  to  station  25  + 00  where  d = 6. 50  (See 
ES-142,  sheet  4. ) Begin  a new  water  surface  profile  computation 
at  station  55  + 00  for  a Q = 1000  cfs  with  an  assumed  depth 
dn_,Q  = 7*97  ft  (where  so  = 0.001).  Compute  the  water  surface 
profile  upstream  to  station  25  + 00  where  d = 6.49  The  differ- 
ence of  the  two  computed  profile  depths  = 6. 50  - 6.49  = 0.01  is 
less  than  the  tolerable  error  of  0.10  ft.  The  last  computed 
profile  is  acceptable. 


Compute  the  pool  elevation  by  adding  the  depth  d and  the  velocity 
'“.21 


head 


v 
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(corresponding  to  the  discharge  Q = 1000  cfs)  to  the 


channel  bottom  elevation  97*00  ft  at  station  25  + 00. 


65+00 


4-16 


That  is, 

pool  elevation  = 6.82  + 97*00  = 103*82  ft  f 

Since  the  computed  pool  elevation  for  the  selected  discharge  is  less 
than  the  given  pool  elevation,  greater  discharges  are  selected  and 
computed  until  the  computed  pool  elevations  bracket  the  known  pool 
elevation. 

Compute  pool  elevations  for  the  discharges  Q,  = 1500  cfs  and 
Q = 2000  cfs,  following  the  same  procedure  as  given  above. 

Plot  the  rating  curve  (pool  elevation  vs.  discharge  Q,)  as  shown  in 
the  following  figure. 


Channel  Discharge  Q,,  cfs. 

Rating  curve  of  pool  at  exit  station  25  + 00 
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Read  the  discharge  Q = 1750  cfs  corresponding  to  the  given  pool  eleva- 
tion. 


r 


Depth  d,  ft . 


Plot  the  rating  curve  (depth  d vs.  discharge  Q)  for  the  channel 
at  station  25  + 00. 
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1000  1500  2000 

Channel  Discharge  Q,  cfs. 

Rating  at  channel  entrance  station  25  + 00 

Read  the  depth  d = 8.65  corresponding  to  the  determined  dis- 
charge Q = 1750  cfs. 


